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ABSTRACT 

FRIIb radio galaxies provide a modified standard yardstick that allows con- 
straints to be placed on global cosmological parameters. This modified standard 
yardstick is analogous to the modified standard candle provided by type la su- 
pernovae. The radio galaxy and supernova methods provide a measure of the 
coordinate distance to high-redshift sources, and the coordinate distance is a 
function of global cosmological parameters. 

A sample of 20 FRIIb radio galaxies with redshifts between zero and two 
are compared with the parent population of 70 radio galaxies to determine the 
coordinate distance to each source. The coordinate distance determinations are 
used to constrain the current mean mass-energy density of quintessence Qq, the 
equation of state of the quintessence w, and the current mean mass-energy density 
of non-relativistic matter Q m ; zero space curvature is assumed. Radio galaxies 
alone indicate that the the universe is currently accelerating in its expansion 
(with 84% confidence); most of the allowed parameter space falls within the 
accelerating universe region on the Q m — w plane. This provides verification 
of the acceleration of the universe indicated by high-redshift supernovae, and 
suggests that neither method is plagued by systematic errors. It is found that 
Q m must be less than about 0.5 and the equation of state w of the quintessence 
must lie between -0.25 and -2.5 at about 90% confidence. 

Fits of the radio galaxy data constrain the model parameter f3, which de- 
scribes a relation between the beam power of the AGN and the total energy 
expelled through large-scale jets. It is shown that the empirically determined 
model parameter is consistent with models in which the outflow results from the 
electromagnetic extraction of rotational energy from the central compact object. 
A specific relation between the strength of the magnetic field near the AGN, 
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and the spin angular momentum per unit mass of the central compact object is 
predicted. 

Subject headings: cosmology: cosmological parameters — cosmology: 
observations — cosmology: theory — cosmology: dark matter — galaxies: 
active 

1. Introduction 

There are several independent ways to determine the global cosmological parameters 
that describe the current state of the universe. It is important to have several complementary 
and independent methods that yield consistent results since any given method could be 
plagued by unknown systematic errors. A particularly useful way to determine global 
cosmological parameters is through measurements of the coordinate distance to high-redshift 
sources. This method is particularly useful because the coordinate distance depends only 
upon global cosmological parameters, and is independent of the clustering properties of the 
mass-energy components that control the expansion rate of the universe, as long as each 
component is homogeneous and isotropic on very large scales. The coordinate distance is 
also independent of whether different components cluster differently, known as biasing, and 
on whether the dark matter is cold, warm, or hot (though it does depend on the equation 
of state of each component). 

Two cosmological tools that are particularly sensitive to the coordinate distance as a 
function of redshift are FRIIb radio galaxies, which provide a modified standard yardstick 
(Daly 1994, 2002a; Guerra 1997; Guerra & Daly 1998 [GD98]; Guerra, Daly, & Wan 2000 
[GDW00]), and type la supernovae, which provide a modified standard candle (e.g. Riess 
et al. 1998, Perlmutter et al. 1999). Some aspects of the methods are compared in §4. 
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The coordinate distance to a source at a given redshift depends on the present value 
of the mean mass-energy density of each component, and on the redshift evolution of the 
mass-energy density of each component. Non-relativistic matter, including baryonic matter 
and clustered dark matter, have a mean mass-energy density that evolves as (1 + z) 3 . 
There is a contribution from radiation and neutrinos left over from the big bang, currently 
negligible, which has an energy-density that evolves as (1 + z) A . There may also be a 
cosmo logical constant, which constant energy density, and so evolves as (1 + z)°. 

The redshift evolution of the mean mass-energy density of a component depends on the 
equation of state w = P/p of the component. There could exist a dynamical vaccum energy 
density called Quintessence (Caldwell, Dave, & Steinhardt 1998). Quintessence would have 
an equation of state w such that the density evolves as (1 + z) n , where n = 3(w + 1) 
(see Turner & White 1997; Bludman & Roos 2001, Wang et al. 2000; or the Appendix 
of this paper). For example, pressure-less dust (i.e., non-relativistic matter) is described 
by P = w = and n = 3; a relativistic fluid is described by mj = 1/3 and n = 4, and a 
cosmological constant is described by w = — 1 and n = 0. 

Two components of the universe are known to exist with certainty: there is a non- 
relativistic component (including baryons and clustered dark matter), with zero-redshift, 
normalized mean mass-energy density Q m , and a primordial relativistic component that 
is negligible for the epochs of interest here (including microwave background photons and 
neutrinos) . 

There must be at least one more component, which remains to be determined. This 
third component could be a cosmological constant, quintessence, space curvature, or 
something else. In the simplest case, there is only one additional component that is 
significant and has a measurable impact over the redshift interval where the coordinate 
distance is used to constrain cosmological parameters and the properties of the unknown 
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component. 

Here, the determination of the coordinate distances to 20 radio galaxies with redshifts 
between zero and two are used to constrain the equation of state and current mean 
mass-energy density of quintessence assuming a spatially flat universe as indicated by 
measurements of the microwave background radiation (de Bernardis et al. 2000, Balbi et 
al. 2000). 

The radio galaxy method is briefly reviewed in §2. Constraints on cosmological 
parameters and the equation of state of quintessence obtained using FRIIb radio galaxies 
as a modified standard yardstick are presented in §3. The supernovae and radio galaxy 
methods and results are compared in detail in §4. The relation between beam power and 
total energy for FRIIb radio jets is discussed in §5. Here it is shown that the underlying 
hypotheses of the radio galaxy method is consistent with current models of large-scale jet 
production in AGN. The results are summarized in §6. 

2. Cosmology FRIIb Radio Galaxies 

The use of FRIIb radio galaxies as a modified standard yardstic to determine global 
cosmological parameters has been presented in detail elsewhere; see Daly (1994, 2002a), 
GD98, and GDW00. Radio images of the 20 sources used here can be found in Leahy, 
Muxlow, & Stephens (1989), Liu, Pooley, & Riley (1992), and GDW00. The radio galaxy 
method is based on the following premises. First, it is assumed that the forward region of 
the radio bridge of an FRIIb source represents a strong shock front, and the equations of 
strong shock physics are applied at this boundary. Second, it is assumed that there is a 
power-law relation between the average beam power of the source Lj and the total lifetime 

a In 

of the outflow t*: oc L- (see §5). Third, it is assumed that, at a given redshift, the 
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sources studied will have an average source size similar to that of the parent population 
of FRIIb sources at the same redshift, where the source size is defined as the separation 
between the radio hot spots. All of these points are discussed in detail in previous work 
(e.g. Daly 2002a,b; GDWOO). 

In addition, in its simplest interpretation, the model works in a straight-forward way 
if the radio power of a given FRIIb radio galaxy does not decrease monotonically and 
measurably with time. There is substantial evidence that this is the case for radio galaxies 
that fall into the FRIIb category, as do all of the radio galaxies considered here (Daly 
2002b). For example, Neeser et al. (1995) find that there is no correlation between radio 
size and radio power. The low frequency radio surface brightness of these sources can be 
explained using adiabatic expansion alone if the radio powers of the sources are roughly 
time-independent (Wellman, Daly, & Wan 1997). Gopal-Krishna, Kulkarni, & Wiita (1996) 
find that a roughly constant radio power followed by an exponential drop in radio power 
would explain the relative sizes of radio galaxies and radio loud quasars in the context of 
the orientation unified model. And, it is shown in §5 that the beam power is independent 
of core-hot spot separation, which would not be expected if the radio power of a given 
source decreases with time. Thus, there is no empirical evidence that the radio power 
of a given FRIIb radio galaxy decreases with time, and substantial evidence that it is 
roughly constant over the lifetime of a given source, followed by a precipitous drop. The 
Mach numbers of FRIIb radio galaxies are about 2.5 to 10 (Wellman, Daly, & Wan 1997), 
implying overpressures of 10 to 200, so the lobe will rapidly expand and drop in radio power 
when the beam power shuts off, as suggested by Neeser et al. (1995). Some authors (e.g., 
Blundell, Rawlings, & Willcot 1999) have studied large heterogeneous samples that include 
sources with all types of radio structure and power and both radio loud quasars and radio 
galaxies. To explain the properties of this heterogenous sample, a model is proposed that 
allows for the decrease in the radio power of a given source as a function of time. There is, 
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however, no evidence that this applies to the select subset of sources (the most powerful 
FRII radio galaxies) studied here, and substantial evidence that the radio powers of FRIIb 
radio galaxies do not monotonically decline over the lifetime of a given source. 

Since FRIIb radio galaxies form a remarkably homogenous population, and all the 
sources at a given redshift are quite similar, the average size of a given source will equal 
the average size of the full population < D > of FRIIb radio galaxies with similar redshifts. 
The dependence of D* and < D > on the coordinate distance, and hence on cosmological 
parameters, is quite different, and 



since < D >oc (a Q r) (see GD98 eq. 6). By construction, the ratio < D > / D* should 
be a constant, independent of redshift. Because the model is based upon a comparison of 
individual source characteristics with the characteristics of the parent population at the 
same redshift, selection effects, such as radio power selection effects, are unlikely to affect 
the constraints placed on cosmological parameters. All of the sources for which values of D* 
are determined belong to the parent population. The sources for which is determined 
are subject to selection effects that are identical to those of the parent population, used 
to determine < D >. Thus, this method is one of the few methods that is unlikely to be 
affected by observational selection effects. Equation (1) allows constraints to be placed on 
the model parameter f3, and the cosmological parameters that enter into the determination 
of the coordinate distance. 



< D > /D* oc (a r) 2 ^ 3+3 / 7 




3. Constraints on Quintessence From Radio Galaxies 



Recent measurements of the cosmic microwave background anisotropy suggest that the 
universe has zero space curvature (de Bernardis et al. 2000, Balbi et al. 2000). Here, a 
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universe with zero space curvature (k = 0) and quintessence is considered. The coordinate 
distance to a source at redshift z in a spatially flat universe is 



where E(z) = ^/fi m (l + z) 3 + (1 — fi m )(l + z) n . The coordinate distance is obtained by 
substituting this into eq. 16 and solving for (a Q r) (see, for example, Turner & White 1997; 
Bludman & Roos 2001; the Appendix of this paper). 

The equation of state w of quintessence is assumed to be time-independent over the 
redshift interval from zero to two (Wang et al. 2000). Thus, the mean mass-energy density 
of quintessence will evolve as (1 + z) n , where w = n/3 — 1. For a present value of the 
normalized mass-energy density of quintessence Qq, and that of non-relativistic matter Q m , 
we have Qq + Vt m — 1 at z — 0, since space curvature is taken to be zero. This implicitly 
assumes that there is only one unidentified component, Qq, that contributes significantly 
over the redshift interval where the coordinate distance will be determined. 

The minimization of the chi-squared of < D > /D* — constant leads to constraints on 
Q m , u>,and the model parameter (5. The one and two-dimensional confidence intervals have 
been determined in the usual way (see, for example, Press, Teukolsky, & Flannery 1992). 

E1 Figure 1 shows the 68%, 90%, and 95% one-dimensional confidence contours in 
the Q m — w plane. At 90% confidence, Q m is less than about 0.5 for any value of w. A 
cosmo logical constant, represented by w = — 1 is clearly consistent with the radio data. This 
Figure differs slightly from the two-dimensional confidence intervals presented in Figure 2 
of Daly & Guerra 2001. At 68 % confidence, Vt m < 0.25, and -1.3 < w < -0.45. At 90 
% confidence, Vt m < 0.5 and —2.6 < w < —0.25. Note, some models for the dark energy 
invoke mechanisms that behave like quintessence with w less than -1 (e.g.Frampton 2002). 




(2) 



NOTE TO EDITOR: Figure 1 should appear here 
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The fifty-four supernovae identified by Perlmutter et al. (1999) and used to obtain 
their "primary fit" C have been similarly fitted, and imply strong contraints on Q m , though 
the supernovae do not place a lower bound on w, as shown in Figure 2. 

E2 

As described by Turner & White (1997) and presented in some detail in the Appendix, 
the expansion rate of the universe is accelerating at the present epoch if 

l + 3w(l-Q m ) < . (3) 

This line is drawn on the Q m — w plane (see Figures 1 and 2); points below the line 
represent solutions for which the universe is currently accelerating, while those above the 
line represent solutions for which the universe is currently decelerating. Radio galaxies 
alone (Fig 1) indicate that the expansion of the universe is accelerating (q < 0); this is an 
84% confidence result. 

E3 It is very important to consider whether significant covariance exists between 
the different parameters determined by the fit. Figure 3 shows the 68%, 95%, and 99% 
one-dimensional confidence contours in the w — (3 plane. Clearly, there is no covariance 
between the equation of state w of quintessence and the model parameter (5: w and (5 are 
independent. 

E4 Figure 4 shows the 68%, 95%, and 99% one-dimensional confidence contours in 
the Vl m — f3 plane. Clearly the two parameters are independent. The data suggest that 
(3 ~ 1.75 ± 0.25 independent of the cosmological parameters Q m , or the equation of 
state w of quintessence. 

E2 NOTE TO EDITOR: Figure 2 should appear here 
E3 NOTE TO EDITOR: Figure 3 should appear here 
E4 NOTE TO EDITOR: Figure 4 should appear here 
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4. Comparison of Supernova and Radio Galaxy Methods 

The observation of the acceleration of the universe through measurements of the 
coordinate distance to high redshift sources is so important that it needs to be independently 
verified using different methods. In §4, it was shown that the radio galaxy method alone 
indicates at about 84% confidence that the universe is accelerating in its expansion at the 
present epoch if the universe has zero space curvature. The supernovae teams report similar 
results at stronger confidence levels (Riess et al. 1998, Perlmutter et al. 1999), as illustrated 
in Figure 2. It is important to consider the similarities and differences between the methods 
to insure that the results of each are independent, and to compare the quantities that are 
actually determined in each method. The two methods are compared in Table 1. 

E5 

The two methods are similar is their dependence on the coordinate distance, and in 
the sense that both provide a modified standard rather than an absolute standard. They 
are somewhat similar in their redshift interval of coverage, and in the number of sources 
studied. The radio galaxy and supernova methods differ in the way the "modification" 
relation is derived; one is empirical and one is based on physical arguments. They also differ 
in that the radio galaxy method is not normalized at zero redshift, and is independent of 
the properties of local sources and of the local distance scale; in fact, nearly identical results 
obtain when the low redshift bin is excluded from the analysis (GDW00). The supernova 
relation is derived at low redshift. 

Thus, the radio galaxy method provides independent verification of the supernova 
results. Current constraints on cosmological parameters determined using the radio galaxy 
method are consistent with those obtained using the supernova method. Clearly, any 

E5 NOTE TO EDITOR: Table 1 should appear here 
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selection effects, or unknown systematic errors, must be completely different for these two 
methods. The supernova method involves short-lived optical events, and the radio galaxy 
method involves much more long-lived radio AGN events. The methods are derived in 
vastly different ways, and are normalized differently. However, similar constraints obtain, 
suggesting each method is reliable. 



4.1. Comparison Allowing for Space Curvature 

It is interesting to compare the constraints placed on cosmological parameters 
assuming that three main terms control the expansion of the universe at the present epoch: 
non-relativistic matter with contribution Q m , a cosmological constant with contribution 
Qa, and space curvature. E6 Figures 5 and 6 show constraints placed by the microwave 
background radiation (Bond et al. 2000), and the high-redshift supernovae project (Riess 
et al. 1998); the results obtained by the supernovae cosmology project (Perlmutter et al. 
1999) are similar and are presented by Daly & Guerra (2001). Each method carves out 
a different part of parameter space because each method covers a different redshift range 
with the supernovae being at the lowest redshift, the radio galaxies at higher redshift, and 
the microwave background radiation at still higher redshift. The three methods are thus 
complementary in their redshift coverage. 

The radio galaxy and supernovae methods are independent of the initial power 
spectrum of density fluctuations, and hence are independent of the index of the primordial 
fluctuation spectrum. These methods are also independent of the Hubble constant, and 
of the baryon fraction. The radio galaxy and supernovae methods depend only on the 
global cosmological parameters Vt m and allowing for space curvature, or fl m and w 



NOTE TO EDITOR: Figures 5 and 6 should appear here 
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allowing for quintessence, where Qq = 1 — Q m , and assuming zero space curvature. They 
are therefore particularly important to study since most other methods depend on several 
other parameters such as the initial spectrum of fluctuations, the Hubble constant, and/or 
the baryon fraction. 

5. Lifetime, Beam Power, and Energy in FRIIb Radio Galaxies 

The application of the radio galaxy cosmology method indicates that (3 ~ 1.75 ± 0.25 
irrespective of cosmological model, as discussed here for quintessence, by GDW00 in 
cosmologies that include space curvature, and summarized by Daly (2002a). The total 
energy released over the lifetime of the FRIIb in the form of a highly collimated outflow 
is E* = Lj t* assuming that Lj is roughly constant over the lifetime of a given source, 
which is supported by lack of correlation between Lj and the hotspot to core distances in 
these sources (see Wan, Daly, & Guerra 2000; Guerra & Daly 2001). Values of Lj and t* 
are obtained independently, using different aspects of the radio data for a given source. 
Two values are determined for each source, one on each side of each source. An indicator 
of possible problems with the model assumptions or selection biases is whether the beam 
powers, total source lifetimes, or total energies depend strongly on the distance of the radio 
hot spot from the location of the AGN. Each of these quantities is plotted as a function of 
hot spot distance from the AGN in Figures 7, 8, and 9. There are 40 data points on each 
figure, two for each of the 20 radio galaxies studied here. 

E7 It is interesting to note that the total energy processed through large-scale jets has 
a rather small range, from ~ few x 10 5 M to ~ few x 10 6 M , though it does increase 
systematically with redshift, as shown in Figure 10. 

E7 NOTE TO EDITOR: Figures 7, 8, and 9 should appear here 



-13- 



E8 

—8/3 

Because U oc L- and E* = Lj t*, the relation between beam power and energy- 
released by an FRIIb is Lj oc £ 3/(3 /3) oc El, where q = 2 to 3 for (3 = 1.5 to 2.0. In terms 
of the total lifetime of the source, t* oc L^ 1 ^ 2 ' ~ 7 / 12 ' ~ 2 /^ for values of (3 of 1.5, 1.75, and 
2.0 respectively. Clearly, (3 remains far from the limiting cases of (3 = or 3. 



5.1. Electromagnetic Energy Extraction from Rotating Holes 

Blandford (1990) gives a summary of the energy and beam power of outflows associated 
with the electromagnetic extraction of the spin energy of a rotating black hole. The beam 
power and total energy available are 

Lj = L EM ~ 10 45 (a/m) 2 B\ M 2 erg s" 1 oc (a/m) 2 B 2 M 2 (4) 

and 

E* = E ~ 5 x 10 61 (a/m) 2 M 8 erg oc (a/m) 2 M , (5) 

for (a/m) « 1, where M is the mass of the black hole, M 8 is the mass in units of 1O 8 M , 
a is the spin angular momentum S per unit mass M: a = S/(Mc), c is the speed of light, 
m is the gravitational radius m = GM/c 2 , B is the magnetic field strength, and B± is the 
magnetic field strength in units of 10 4 G (Blandford 1990). 

Empirically, E* oc ^ 3 with (3 = 1.75 ± 0.25. This is consistent with equations (4) 
and (5) when the magnetic field strength satisfies 

B oc M^-V' 2 ^ (a/mf'^-V . (6) 
This is particularly simple when (3 = 1.5; in this case B oc (a/m). For the cases (3 = 1.75 and 
E8 NOTE TO EDITOR: Figure 10 should appear here 
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2, the empirical constraint is consistent with equations (4) and (5) when B oc (a/m) lA M 0,2 
and B oc (a/m) 2 M 1 / 2 respectively. 

Typical values for the beam power, total energy processed by a jet, and lifetime are 
Lj ~ 10 45 erg s _1 , i? t ~5x 10 5 M o c 2 , and t* ~ 10 7 yr ( see Figures 7, 8, and 9). Equations 
4 and 5 imply a total source lifetime of t* ~ 10 9 /(M 8 _B 2 ), indicating a lifetime of about 10 7 
yr for M 8 Bj ~ 10 2 . This is satisfied for M 8 ~ 10, and B A ~ 3. In this case, the beam power 
is ~ 10 45 erg s _1 for (a/m) ~ (1/30), and the total energy is ~ 5 x 1O 5 M . These values for 
M 8 , £? 4 , and (a/m) seem quite reasonable. The scaling between variables required by the 
empirically determine relation between total energy and beam power are discussed above. 



6. Conclusion 

Type lib radio galaxies and type la supernovae are particularly important methods to 
develop to constrain the global cosmo logical parameters Q m , and space curvature, or 
fl m , Qq, and the equation of state of quintessence w. These methods depend only upon 
the properties of global cosmological parameters, and are independent of other factors such 
as the index of the primordial power spectrum, the Hubble constant, the baryon fraction, 
the properties of the dark matter that clusters around galaxies and clusters of galaxies, any 
biasing of dark relative to luminous matter, etc. The radio galaxy and supernova methods 
are completely independent and have completely different potential systematic errors, as 
discussed in §4. 

Radio galaxies may be used to constrain the mass-energy density and equation of state 
of quintessence. These results are presented here assuming a spatially flat universe, which 
is supported by recent measurements of fluctuations of the cosmic microwave background. 
Radio galaxies alone suggest that the universe is accelerating in its expansion at present 
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(see S 3), consistent with results obtained by the supernovae teams, indicating that both 
methods are working well and probably are not plagued by unknown systematic errors. 

The implications for models of energy extraction for the cases of an outflow related to 
accretion and the Eddington luminosity, and electromagnetic energy extraction of rotational 
energy are considered. If the outflow is produced by the electromagnetic extraction of 
energy from a rotating black hole, then the magnetic field strength must be related to the 
spin angular momentum of the rotating black hole S, the mass of the black hole M, and 
the gravitational radius of the black hole to, as described in §5. The relation is particularly 
simple if j3 = 1.5, and implies that the magnetic field strength satisfies B oc (a/m), where 
a = S/(Mc). 
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our studies of radio galaxies. We also thank Megan Donahue, Paddy Leahy, Chris O'Dea, 
Adam Riess, Max Tegmark, and an anonymous referee for helpful comments. This work 
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National Science Foundation, and by Penn State University. Research at Rowan University 
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Table 1. Comparison Between Supernova and Radio Galaxy Methods. 



Super novae 



Radio Galaxies 



Type SNIa 
oc {a r) 2 - 
< z < 1 
~ 100 sources 
modified standard candle 
light curve =>- peak luminosity 

empirical relation 
written in terms of observables 
model parameter determined at z = 
may depend on local distance scale 
universe is accelerating 

some theoretical understanding 
well tested empirically 



Type FRIIb 

oc (a r) L6 
< z < 2 

20 sources (70 in parent pop.) 
modified standard yardstick 
radio bridge =>- average length 
physical relation 
written in terms of physical variables 
model parameter not determined at z = 
independent of local distance scale 

Q m is low 

universe is acclelerating if k=0 
good theoretical understanding 
needs more empirical testing 



-19- 



Fig. 1. — One-dimensional confidence contours (68%, 90%, and 95%) on the Vt m — w plane 
obtained for a spatially flat universe. The dark dotted curve separates the plane into an 
"accelerating" and a "decelerating" region. The best fitting values of parameters, as well 
as the chi-squared of the fit (obtained for 20 radio galaxies and 16 degrees of freedom) are 
indicated on the figure. 

Fig. 2. — One-dimensional confidence contours (68%, 90%, and 95%) on the fl m — w plane 
obtained for the 54 supernovae in the "primary fit" C of Perlmutter et al. (1999) for a 
spatially flat universe. The dark dotted curve separates the plane into an "accelerating" and 
a "decelerating" region. The best fitting values of parameters, as well as the chi-squared of 
the fit (obtained for 54 supernovae and 50 degrees of freedom) are indicated on the figure. 

Fig. 3. — One-dimensional confidence contours (68%, 90%, and 95%) on the w — (3 plane 
obtained allowing for quintessence. Figure shows that (5 and w are independent. 

Fig. 4. — One-dimensional confidence contours (68%, 90%, and 95%) on the Vl m — (3 plane, 
obtained allowing for quintessence. 

Fig. 5. — The two-dimensional confidence intervals in the Q m -QA plane from the cosmic 
microwave background and FRIIb radio galaxies. Solid lines indicate the constraints from the 
cosmic microwave background (Bond et al. 2000) and shaded regions indicate the constraints 
from FRIIb radio galaxies (GDW00). The combined constraints indicate that there must be 
a significant contribution from a component other than non-relativistic matter at the present 
epoch. 

Fig. 6. — The two-dimensional confidence intervals in the Q m -Q\ plane from supernova and 
FRIIb radio galaxies. Solid lines indicate the constraints from the high-redshift supernovae 
team (Riess et al. 1998) and shaded regions indicate the constraints from FRIIb radio 
galaxies (GDW00). The 1 a concordance region is intersected by the k=0 line, suggesting 
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that the universe is spatially flat with a low value for Q m . The results obtained by the 
supernovae cosmology team (Perlmutter et al. 1999) are very similar to the supernovae 
results shown here. 

Fig. 7. — The beam power, obtained by applying strong shock conditions to the forward 
region of the radio source, is plotted as a function of hot spot distance from the origin of 
the source. Two values are obtained for most sources, one for each side of each source. The 
normalization is described in WDG00 and GDW00. Filled symbols are the six new radio 
galaxies obtained from the VLA archive and described in GDW00; open symbols are the 14 
original radio galaxies described in WDG00 and elsewhere. 

Fig. 8. — The total time a particular AGN will produce each of its two large-scale jets is 
plotted as a function of distance from the origin of the source. The filled and open symbols 
are the same as in figure 7. 

Fig. 9. — The total energy a given AGN will process through each of its two large-scale jets 
is plotted as a function of distance from the origin of the source. The filled and open symbols 
are the same as in figure 7. 

Fig. 10. — The total energy a given AGN will process through each of its two large-scale jets 
is plotted as a function of 1+z. The filled and open symbols are the same as in figure 7. 

7. Appendix 

The deceleration parameter is defined to be q Q = — (aa)/d 2 evaluated at z=0, where 
a is the cosmic scale factor. This may also be written q = —a /(a H%), where quantities 
evaluated at z=0 have a subscript 'o', and H Q = (a /a ). It is straight-forward to show that 
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(e.g. Peebles [1993] eq. 5.15) 

(a/a) = — ^(p, + 3 Pi) = 3- X^ 1 + ( ? ) 

where pj is the pressure, p, is the mean mass-energy density, «)j is the equation of state of 
ith component, Wi = Pi/ Pi, and the quintessence or cosmological constant term is included 
in the summation. 

For each component, mass-energy conservation implies 

Pi = -3(pi +pi)(a/a) (8) 

(e.g. Peebles [1993] eqs. 5.16). By definition, the equation of state is Wi = Pi/pi, so the eq. 
8 implies (p\/ Pi) = —3(1 + Wi)(a/a) . When the equation of state is time independent, 
the solution to this equation is pi = Pi iD (l + Z ) 3< - 1+Wi \ where (1 + z) = a Q /a. Thus, a 
component with equation of state and present mean mass-energy density p Q will have a 
mean mass-energy density at redshift z of p = p Q (l + z) nt , where rii = 3(1 + Wi). 

When k — 0, 

(«/«) 2 = yE^ (9) 

where p t = p i:0 (l + z) n . It follows that Hi = (a /a ) 2 = ( §2 f-)J2Pi,o- For k — 0, Pi,o = Pc,o 
where p Cj0 is the critical density at redshift zero. Thus, 

Hi = -3- Pc,o (10) 

for k — 0. The deceleration parameter q Q = —a /(a Hl) then becomes 

q = (1/2)^^(1 + 3^) , (11) 

where fij = Pi, /p C)0 and equations 7 and 10 have been used. 

The universe is accelerating rather than decelerating when q Q < 0. This can only occur 
if (1 + 3wi) < 0, or Wi < —1/3, which is a necessary but not sufficient condition to have 
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an accelerating universe at the present epoch. If there are only two significant types of 
mass-energy controlling the expansion rate of the universe at the present epoch, quintessence 
and non-relativistic matter, then the deceleration parameter is q Q = VL m /2 + £Iq(1 + 3w)/2. 
The universe will be accelerating in its expansion when 

l + 3w(l-Q m ) <0 , (12) 

which follows since VLq = 1 — Q m . Thus, a curve can be drawn on the Q m — w plane; points 
below the curve represent solutions for which the universe is currently accelerating, while 
those above the curve represent solutions for which the universe is currently decelerating 
(e.g. see Figure 1). For example, this curve is shown on Figure 1, which illustrates the 
constraints obtained using radio galaxies alone, and is discussed in more detail below. Note 
that radio galaxies alone place interesting constraints on Q m and w, and these results are 
consistent with those obtained using other methods (e.g Wang et al. 2000). 

The coordinate distance to a source at redshift z follows from the equation 

J dr/Vl - kr 2 = J dt/a(t) = (l/a a ) J (a/a^dz (13) 

(Weinberg 1972). For a spatially flat universe, the left hand side of the equation reduces to 
the coordinate distance r. Equations 9 and 10 imply that (a/a) 2 = (H 2 /p C:0 ) Pi :0 (l + z) TH 
or 

(d/a) 2 = H 2 J2m + zr ■ (14) 
Following the notation of Peebles (1993), this reads (a/a) = H Q E(z), where 

E(z) = y/^m + z)^ ■ (15) 

Thus, the coordinate distance to a source at redshift z in a spatially flat universe is (see 
equation 13) 

(a r)=H- 1 f dz/E(z) . (16) 
Jo 



-23- 



A spatially flat universe has = 1- The components that must be included 

in equation 15 are those that contribute from redshift zero out to the redshift that 
the coordinate distance is being determined. Two components are considered here: 
non-relativistic matter with normalized mean mass-energy density at z = of Q m , 
and quintessence with normalized mean mass-energy density at z = of Qq. Thus, 
fl m + Vlq = 1, and equation 15 becomes E(z) = y/VL m {l + z) 3 + (1 — fi m )(l + z) n . The 
coordinate distance is obtained by substituting this into eq. 16 and solving for (a Q r). 
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